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FOREWORD 


Tnis is the final report submitted in fulfillment of Contract 
F04611~-67-C-0019. The program had two Air Force Project Engineers, 
Lt. James 0. Cox/RPRRC and, later, Capt. J. F. Ensminger/RPRRC. 

Volume I of this report contains the background, development, and 
critique of the Sensitive Time Lag Theory of combustion instability, 

a description of the computer program that is used to provide a solu- 
tion of the analyticul model, and the use of the model in the analysis 
and design of liquid rocket engines based on empirical] correlations. 
Volume II contains the flowcharts and listings of the computer program. 
Therefore, Volume I is referred to as the Engineer's Handbook while 
Volume II is referred to as the Programmer's Handbcok. 


The contract effert was conducted by the Thrust Chamber Engineering 
Division of Liquid Rocket Operations under Dr. C. B. McGough, Program 
Manager; Dr. R. J. Hefner and, during the final phases of the contract, 
by Mr. J. M. McBride, Project Managers; Mr. A. J. Smith, Jr., Project 
Engineer; Dr. F. H. Reardon, of Sacramento State College, and 
Dr. L. M. Crocco, Dr. W. A. Sirignano, and Mr. D. T. Harrje, of Princeton 
University, Consultants. This report has been prepared in accordance 
with MIL-STD-847 (USAF) dated 25 February 1965. The period of this 
report is 2 September 1966 through 31 November 1967. 


It is essential to give special acknowledgement to Dr. L. M. Crocco, 
originator of the Sensitive Time Lag Theory, and his co-workers whose 
individual contributions to this discertation are reflected throughout. 
In addition to contributing significant portions of the Background, 
Theory, and Conclusion sections, as well as Appendices I through III, 
they have reviewed the entire text and have made many helpful suggestions 
during the entire course of the contract. 


This technical report has been reviewed and is approved. 


Capt. J. F. Ensminger 
USAF Project Engineer 
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ABSTRACT 


The main objective of this report is to include, under one cover, 
all of the work concerned with the development of the Sensitive Time 
Lag Theory of liquid rocket engine combustion instability. This work 
tncludes all tiie aspiztc 7* the theory from the mathematical formulation 
of the analytical modei to the application of the model to actual engine 
problems. The initial section of the report reviews the logical cen- 
siderations of the instability phenomenon and relates how the time lag 
concept conforms analytically as well as experimentally to the problem. 
Thereafter, the mathematics of the model are developed with the major 
emphasis placed on the linearized model; however, various aspects of 
the nonlinear model are also discussed. The mathematical analysis gives 
rise to a computer program, which is presented in Volume I in the form 
of an Engineer's instruction manual and in Volume iI in the form of a 
detailed description for the Programmer. The report then focuses its 
attention on the designer and instructs him in not only how to use the 
model in practical situations but also how to interpret and correlate 
test data. The main body of the text concludes with a critique of the 
time lag concept and outlines the kind of research that is needed in 
order to imprcve the time lag model. 
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mY @eO1, 
SRTOQSZIO 
INTAOPI0 
TRTERVC® 
1NTEaOS0 
1m 100940 
InTeed7O 
IHreoveo 
14TSO990 
TFRT619000 
INTO1610 
14761020 
imy61es¢ 
INT61640 
1RT616$0 
THYS1060 
THTE1070 
1761000 
IMTO1G699 
TRT61190 
IMTS11I0 
IMT@L220 
TR161180 
TmT@}, 40 
Twre.p er 
14521100 


BREA ORSATE SOLIS" MEF DEC E oM/90 MELMOD 
° 
. Date -- TODav-* DATE Im Te FORM OO mmeeH vie 
* “Are TATE (MCDA Lhe WHERE BCOATE £5 2 CELLS WHICH WILL COMTatl 
rz “mE BCD FORM OF TODAY-S Oa E 
. 

Ewtay cave 

TSx ENTRY 


oon 


OECe 
2e mee, 
le yates 


NCre 
FEB. 


SIOFTC CNTs LIST Moe 


OCT. 


SUAROUTINE INTSIX ev oKT 9D} 
DIMENSION KCODOV COTO RCLADVC (A! 


EQUIVALENCE 


SUILO RETURN ADORESS 
TSU ERTRY 


«0aD LP XRe FOR CALL SEC. 


—™m00YY 
mM /ODYVOO 


Commons YYCOOC 
D00006/YYD000 
000000 ‘00907 
060 vY/0000C0 
00197V/000000 


0/BAAAsAB0000 
DOCOOOSBAAA IM 35a CIMIIRECTS 


480000/000006 
ABOO0Oe YY 


ESCAPE 


SEPT AUGe JULY JUNE MAY 


APR, 


echa 
FDA 10 
*eDa 20 
erOA 36 
onda 40 
ero, 89 
Sx5a $9 
cos 76 
oA 8D 
“roa 90 
*LOA 
@KOA 
nA 
*nOa 
KDA 
onca 
enna 
eK OA 
omy 
*EDA 
KOA 
*KOA 
noe 
*KCA 
POA 
enpa 
ecOa 
@KOA 
ecoa 
eKOA 
KDA 
#KOA 
*KDA 
KDA 
econ 
eKOA 
*KOA 
KOA 
exOA 
exOAa 
@LOA 
KDA 
xda 
eKOA 
PKA 
®KDA 
*KOA 
KOA 
ROA 
ex0a 
@KDA 


INTs 
INTa 10 
HTS 29 


EXCEL DMD AIKCLZ[ eX oCKCCIF oMIFOUKC(MIORSPSLVCULISYRIINTS 39 


Tel VEE eVIF EMME CAPEK et VCs eVa) tNTa 45 
19 ASSIGM 9% TO NB Inte 36 
Je2 IRTe 60 
Sext INTRA 76 
20 1F1X039990 260030 IMT® 8&0 
30 GO TO MA2( 900389) INTe£ 96 
40 TFCYE 31130050030 TINTS 100 
$0 1F 1 I-2160060+70 tNT& 110 
60 YEG, TINTS 120 
GO TO 180 tHTe 130 
7) ASSIGN $30 TO WA INTO 140 
Je le} IN,& 75C 
BO Klexi yy) TNT& 160 
K2exX( Jad? INTS 170 
RYeon( y-2) INT& 360 
vievey) TINTS 190 
Y2evtl=b) INTS 200 
Viev( ja?) INT& 210 
GO TY WB011300170) TATe® 225 
OO TFEK J) "8112001999190 INTS 230 
100 TF 322113001300110 TINTS 240 
210 ASIGN 140 TO Na INTe 750 
1270 Jase INTS& 260 
4H TO 20 INTS 270 
130 09 HM Japed 1%Te@ 260 
Ke syakO 2) INTS 290 
WeO VCE SDeV INTS 300 
150 DeK2-x} tNT&e 310 
Ay eR-X] INT 320 
A2*8-X2 INT’ 330 
VECALOAZ/Z eO/DOLLY3—VIV/ERIRANZ HCV IVI I/D RAZ/D BV IGAL/0EV2 INTS 340 
GO TO 180 INT& 350 
180 ASSIGN 170 TO MA INTA 300 
GO TO 80 INTe 370 
170 Keexty-4) INT&S 380 
Yeey(Jj-3) INTS 390 
Devrnax? INTS 400 
Ay eA-X2 INT& 410 
AP eR-X3 INTS &20 
MMYQ2CCVInVLP/IKI=x1) INT® @30 
KM2301Y3-¥29/0 INTS 440 
XMYROTYE“VVI/(KG-X3) TINTS 450 
VER) PA2FEI/2 O/NSS2OTAK) 2-KM2 9) OAZOA1EO2/ 2 oO/DSO2E EMIH—KM23}V-AZPINTE B60 
1V2/D4al#V3/0 InNTa &70 
180 YOuove TINTS 880 
RETURN INT& 490 
EAD INTs $00 
STAR IC MTane {1ST Moe TINTON 
SUMBOITINEG INTAOT Xe Vent oVOsT? TNTan 10 
QUMENSTOM MEDb eT EOb AKC (ATSC fa? TNTAN 20 
EOITVALEMCE EXCEL KEP OE MOLD oP eCKCE VIM SKCH oC eCKCLL I SYIVINTAOD 30 
DsEMC EDD AVDA ELVOE IAB CKO at rat INTAaD 40 
10 ASSIGN 99 TO AA INTAD $0 
Jo? INT&™ 80 
Aaxt INTA&O 70 
oO EFIMTIIIIOS6O030 INTED 80 
30 GO TH) Nhe! 200180" INTOD GO 
AD TREE FEDS OF OIN THT&5100 
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fees $reg 


ent er: 
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30 
60 


70 


60 


9 
100 
1320 
170 
136 


146 
190 


166 
270 


180 


TF ( 2166060079 
Y¥e0.0 

GO TO 180 

ASSIG® 155 TO Ae 

Jel 

x1exK(s) 

X2eX(J-13 

K3eK(Jo2) 

Yyevts) 

Y2ev{ I-31! 

¥3eVt 5-2) 

GO 15 Oe'1S05170) 

TF (X1I9-B122992000100 
TFC 5-2929002900210 
aSSTGR 166 TO Wa 

jeje} 

© TO 20 

BO 140 Jeles 
CHIEN) 

vet sre¥tar 

Dex2ni 

PleB-A} 

AZ=B-K2 
AMIPOCVI$V2)/0K9-¥2) 
RMLZELVI-YLIZURT=RI 
RM2H0(KN29-E129/200/0 
VOCALOAPORMIB-AZOY] /OOA19Y2/0 
DVeRMIBO (A; OA? 1¢xMIZ 
GO T0 180 

ASSIGM 170 TO AB 

e& To 80 

Keoki J-3} 

YaoY(J-+3} 

Dexge-X2 

AleB-X2 

A2eB-K3 

MALZOUY2VVP SIKH) 
LMPBOMVS-¥20/0 
RMyReTYRHV3I/(Xa-x3} 
AMZOATO(KM) 2—-EN23} 
AALOAL OL ENSHeKMT9) 
YO~(A1992/200/D8 1 AMZ OAM) } -ADOVZIOA16791/0 
OVet AMZO(2.OOA10A2 }OAM10( 2 oGOA20AL) 1 /200/DOC2 OEMS 
@eTURN 

£40 
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40 


$0 
66 
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86 


SUBROUTINE PAGEILIMES! 
WEAD MOVED TO /PROLOG/ AND PAGE MODIFIED 10 PRINT HEAD 255 47 


COMEON /PROLOG/ LOGIK(SSSs MEADILZ2}0 Side Size EORI 
OIMENSIO"N TODAY (2) 
OATA EP6 / 0 7s 


TF UL IMES=60520020010 
Lez 

60 10 & 

KoLon tags 
1F{K~60190020990 

tek 

RETURN 


LeLInesez 

TF ¢ KPO EO. 07 CALL DATE (TODAYS 

KPG = KPG e 

WRITE (60801 TODAYs HEADe EPG 

FORMAT { IH] Si GMDATE = =2A60 12 12A6e 31K SHPABE 15 1 
@0 70 40 

tno 


IMTAO11C 
THT40120 
IMT80130 
IRTOOE8O 
1mT40180 
IMTOO 60 
1MTAO170 
ImTaO180 
[mT4D190 
1mTad200 
INTAOZI0 
INT40220 
INTAO2306 
ImTaD249 
Imveo2S¢ 
tmTaQ0260 
IMTaD2Z70 
TATAO2E0 
INTADI90 
1nNTs0300 
1NTa0310 
IRTaV$20 
INTAOS30 
INTao3eo 
ImTaosse 
INTeO960 
INTAD370 
TnTao3so 
TaTeD990 
1aTeDeee 
INTapele 
taTaoare 
» 700030 
140180000 
1NTA0a50 
InYe@ega 
InTepere 
1nTe0086 
INTADeto 
1nTe0500 
taTeos10 
1NTA05208 
TRTZSBIC 
TRTOLI00 
INTADS 40 


SOP IGIN SLO0E 

SINCLUDE aaror 

SIBV-C CHa LIsTomes Cham 
SUSPOUTINE CHAMBR Cua =610 

i cua = 20 

¢ 20 $5P &Y MODIFIED FOR TASULAR INJECTOR COEFFICIENTS CHAN 30 

€ cuam 40 
LOGTCOL LOGE Ce ARUN BRN oC RUN DRUM o£ PUN CFRUNSGRUNCHRUN. I RUN. RUN CHAM SC 
LOGICAL SLie *L2¢ EOS Cham 66 
Ata, “acH cua FC 
CCewe 4s cHam 680 
1 Game wate WITEIODs AVNTIOe BYNCIOTs CVMRCIOD> CVMELIONe EESDSCCMAM 90 
COMMON §«/PROELOG/ §LOGIKRI38}s HEADI}2)¢ Site SL2e EORY cuad 100 
COMBEDN/ABCOF / 01 NP 6S TCe cat 320 
DIMEMSTOM EXTRAI}COle welts) cyan 120 
DIMENSION OIG LO3OO) SALTS oBID PCs sd eDET? Cran 330 
DIMEGSIO—m «KE 1397 o V8 1335 006294) eSTOWL 222) oSTODATI 06071 CHAM 340 
Dimensior 221203) Cuan 130 
OEMENS IDG GIL) DISTLE207 OSSTMEZO} Cham 369 
Clegastin asit¢sor cuant 170 
EQUIVALENCE { EXTRAs Simms STODAT 1 CHAM 140 
Ou! vaLENce CHAM 150 


1 {LOGIECL) © ARUN: (LOGIRC2) © BRUNTe ff. 0GIE69) » C@umte CHAM 20C 
2 (LOGTREO) © ORUNIs  CLOGIKRISS 6 ERUMDe (LOGIKIO) « FRUMI. = «CHAM 210 
%  CLOGIRCTD © GRUMD> «FLOGIKIGD » HRUMEs  CLOGIKI9) . EMUMie CHa 220 
ry 


{LOGIEI107¢ Jt? CHa 230 

EQULY LEMCE SEXTAGES Do AD eCERTAALZ2) oCOrelENTAALSIOCCHSCERTRAIAIs CHAM 200 

1CD) oe CERTRAIS) cCEDe EXTRALG) oCF 1 o CEXTRALT 1 0CG) o(EATRALS) CHI cman 250 
EQUIVALENCE (EXTRAG9) 0CId Cham 260 
EQUIVALENCE (D1MP613 0A) 2 (DIMP (S003 } Ble (DIMPIS001) Die Cham 270 
1LDSP E3903 2 G2 c(DIMPCSBOLE CI oMERTRACZID CWC) © (BI ebeUE) CHa 280 
EQUIVALENCE (O10P 1601 1c AMET) c (D1NP 1281918 00% co (DINPL 39011 eKde Cuan 290 

! (Osh 0545) Devt CMa 300 
€Oulvacenct Chas 310 

1 CSTOWES DTZ Ie (STOWE 2IZboVMD CHAM 320 

2 6 ESVOMEZIS Toth e CSTOMLZISD MED e (STOMEZIOD VLD (STOWEZITD KIICHAN 330 
Do «*STCWE2EED MDs CSTOWIZI9D oKERD + (STOWL220) oAN) cHan 340 

& © = €8T0m1 27229 KQUADI > *STOWIZZZIOKLD Cuan 350 
So «= KERTRAUSE De OLSTLD© CERTRALTIN» OFStM D CHAN 360 

& 6 = KEXTRATLON Cade LERTRACZ}1 oSMM te CEXTRACIZ} > MACH? Cuan 370 

7 ° (03 P 1207) U1 BAR) CHAM 380 

C Cuan 390 
C9900090000660 000005006868 oO HDOFERSOES OOF SE CODES OOECEOCOEES OD ODLOER COED CHAM £00 
¢ coam 610 
10 FORMAT TINO 60H a 8 c ® € € G " CHAN 420 
1t J = Ff 3WotlOFG.01) cwan 430 


20 FORMATI//7 69K 0 1OQHSCESSeeSLeZe—RseCeTeSOOS THE FOLLOWING MAIN CONTCHAN 640 
IROL DATA WILL BE USED 1% THIS CASE svesveccaooooooooosocese § 64///sCHAM 850 
ZASKoDSMMATIO CF SPECIFIC MEAT 1GADMA) © oF Tobe//s4SKo22HOESIRAIO MACHAM SOC 
SCM HUMBER © 6612650266 (20 SF BEING CALCLLATEO) £/49% LTecHCHAN 470 
SAMBER RADIUS © oF Te3e OH CINCHESI 0/7/0450) THCHAMBER LENGTH © oF ToXCMAN 480 
be OH LIMCHES) 77 045K+ LIMSPEED GF SOUND © oF 10030 Ot IFT SSECD0//0CHAM 490 
TOSXe ZIMCHAMBER MODE DESCRIPTION © of Se5e 20H 100 FOR LONGITUDINALCHAM 500 


B MODESH eso) cHam $10 
30 FORMATI///o33N pbBHERCCCCOCES CHAMBER FREQUENCIES (WC) eeeceecons Cram 520 
terse) cuam $30 
40 FORMATISF20.$e/ } Cham 940 


50 FOMMATISONs GAHOBCOD MACH AISTRIBUTION IM CHAMBER AS A YULACTION OF CHAM 350 

1 LENGTH ©8066 6// 013 Ko THCHAMBER o 1 SK oGHMACHe 14% e THCHAMBER 1S Xe tHMACMCHAN 560 

2914 e THCHAMBER 915K oSHMACH o/ 0 L2XoGHLENGTMe LIXo1ZHDISTRIBUTIOMSIIRe CHA STO 
DOMLEMGTHs LIX oD 2HDISTRIBUT 1OMel IX eGMLEMGTHS 1 IXOL2MDISTRIGUTIONS/¢) CHAM 580 

60 FOMMATIOL1ON oF 300511 cwam $90 
€000000000000000055 00000000000 uCeTOOOe CESEEOREHOSESOCSOPCOSLCOENOROOES® CHAM 600 
COCOOOOLODOFEDODOSEOELELE SORE LEOSELOTECOOLESEEESENERSOCRDESCCONLECRCES CHAM 610 


id cna 6290 
Fs S$e21 GO TO 175 CHAM 630 

Cc CHAN 660 
$L2 © o TRUE. CHAM 630 

00 70 I © ts 4300 CHAM 660 

70 OLN £1) © 0.6 cHan 670 
GO TO 110 CHAM 680 

¢ CHAM 690 
C0000080065000003005000000000500 000 CSc SESS TOODECOSCESEECEETESECCCSEECe CHAN 700 
ig cwam 71h 
60 WRITE (69100) WE CMaM 72> 
90 CALL EXIT CHAM 730 
190 FORMAT (IMOLOXITMINPUT ERROR: XE = [3 19M, WEMCE TERMINATION } CKal 760 
cMam 750 

09000000 00000050000 10000 00000000 0000000080 00CSSSETERELOTONSSSESES CHAM 760 

oe: SFRSCROSOHOSOOSSSSSSSSEOOHSSSHUSICVOSCOSSS SSsesecoesesesaase CHAM 770 

CHAN 780 

110 KER © 0 cham 790 
CALL OVCHR CK CHK) cwam 600 

¢ CHAM 810 
00 1:6 Te le 10 CHAM 829 

120 DIaPttT) © 0.0 cua 830 
CALL ASTIO ¢ DINPL1)9 HEADL LP» NE D CHAM 64¢ 

IF ( ME oME. 1:9 GO TO 60 CHAM @50 

ARUN ¢ CA 000 cham 860 

BRUN © CB 0.0 CHAM 876 

CRUM © CC eNEe 000 CHAM 660 

ORUM © CD eNEs 000 cMam 890 

ERUN © CE ONE 060 CHAM 900 

FRUM © CF UNF, O6f cham 910 

TRUM © CT oNFe M60 CHAM 920 

IMUM = CJ oMEs 060 CHAM 930 
OGRUN « OIMPT22) oLEe O20 CHa 040 
eANDe 1 ARUN sORe ARUN ORs CRUNM cOR, TRUM) cham 950 

CHAM 960 

MYT MEW MAIN CONTROL CHAM 970 
soeeonecesee @eeseoseeees CHAM 980 

CALL PAGE ¢ 601 CHAM 990 
WRITE (6010) (DIMPETDs Te1010) CHAM}OCO 


WRETE (60291 DIMPIT11s MACHs DINPI 1610 OINP (1530 DINPIi6)e SNH CHAMIOIO 
COSRSOSEESSOSE SCOTS EOTECEOHOTESSH DASHES SSESOSETHTESEROSCEESFUSCEREOCHESOS CHAMIO2Z0 
c ARE FREQUENCIES TO BE CALCULATED CALCULATE AND PRINT. CHAM1030 
COOK ROSSOCR OOH OSEREEDES ORCS O DES EESESEROSESSOS CHAMIOSO 


IF (¢ oMOT. GRUN } 40 710 130 CHAMIOSO 

CALL CEMMEG ¢ WOTL) CHAM1060 

IF « BRUM .OR, ARUN } CRUN © oTRUES CHAMIO70 

1f ¢ CRUN ) ce # CC + 1100 CHaM1 080 

130 JOMEGASABSTEKTMA122014%22,0001 CHAMI090 
WRITE (030) CHaMy100 
WRITE (OeGOVTENTRALIT of 023+ JOMEGA) CHaAM}110 
WRITE to050) CHAM}120 
WRITE (OoHO1f DISTLE THe OISTMITds DISTLITeThe O1ST4¢567)CHaM1130 

Le DISTLEFetade OLSTMtToetaty Feiee te MISTLITI. DIStTMI7)s CHAMIIEG 
POISTL EDO) s OFSTMlyay CHAMI150 
CESESOOOESSSSSFSSESOSOLEEE HED /EOTESUSSESECSSEESESVESESSEOCESSSSESRESE SEED CHam1160 
¢ CHAMI}70 
EOR} © FANN WOR. JRUN .ORe [RUN CHAM1180 
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IF ¢ Cum OR. EORy ! CML PAME ft 40) 
100 1F ¢ MOT. FORD } 10 200 
WRITE (103 stooat 
SACKSPACE 18 
ct « CE 
ost © CI 
Gan © OTM I11> 
10 ¢ .ROT, 18m 1 @ TO 160 
mule t JOMEGA22 1/2 <3 
00 190 fo 3. aul 
190 UITITd © OLMPE 205e2) 2 
wiTecat) © O1NPC JOMEGA } 


¢ 
166 Attu 

c eoeeee 

c PROSANM RETURNS 1 AND ONLY 17 INJECTOR ORBGRAN 15 10 GE RUN 

c 


370 READ 234) S10DaT 
BACESPACE 14 
Ff ¢ mot, Camm } @ 10 200 


0587694083 © Met 
00 100 Te je mul 
O31 104922) © WITtT) 
BIMPL 1003998) © ave 11) 
OmP 1100996) © OVE II) 
¢5eaSTSs © Cymetis) 
160 6900900) © Cuntits 
30 1 1 © We Cee 17 
190 DIMM 13009239 © 0.0 


¢ 
COPPLOLG SEOSBOLGESOG SOLO CLLO COTO IAEOS OOS CO OODIOLISI OHOGOSHPMEIELIOOLE 
¢c 
200 Wiz) © ABSE UCI2) D 
ad 


COQSPSLCOSSSESES HT SESSDOEOCSHLSTOOOLSS SS00FTL OS CEPT CHO MESES ET ESTES: 
c FUN MBER OF FREQUENCIES AMD TEST OESIME POR iQTEMMAL C FLAG 
CETOLOC PSU SEDSAESPECTOS OSE 7900069909 05 00 CHE 50 0SDISISSOSESE COST IETS ONE 
220 MONES—UC 12 1¢.000) 
I CENTAALL2) 129002200230 
220 0107 190021¢3.0 
60 To 200 


200 17 « .meT. “mm @ Te we 
(086080 099000900 20880 0820000006 5000028 BUOIOTLITIETS FORT DOSS POTROOO OTTO 
i¢ SET WH DATA FOR PROGRAM Ce. 
(9990600066000.060008 08880000000 0604000000890 00S OT0TISSTOSOS0STOEOS 9000 
250 D1 (3601 po ERTAACI15 
wee © ROMEG/2 « i 
2241 :e@xTRAL 14) . 
2U(ZICENTRALIS) 
2USSCERTRAL IG) 
O19 (9003) © EXTRAII4) 
5113803) semar 
lg oe e€8. 900 } O1MP(S8EF) = 16108 
e 
ZAVE*( (ERTRAC11 163 001/2cODOSERT {OF GP 9000 1 CRIMP 1 908711 
SP CEXTRACZ11 927002600278 
260 $M0Z © 0.0 
CEOOOO9EESE-0F0-08-00 0808 L 1000090000080 FEELS C8 SCOODESCOCOSISOSOSESSCSCLOSEO 
¢c USING MALF OF CHAMBER FREQUEMCIESs CALCULATE NOZZLE PREQUENCIES 
€ FOR USE Im PROG C. 
COOBOSS 009 OEOOESSSSF SO SSL OE SHEEHDOSGS COS OSOLOESOOOOSOS OSE OS OEOE OSE SOOO 
Rom © EXTRA(15» 
60 tO 260 
270 ROM © EXTRAL1O} 
GRAD © FORTE (2 00/ (EXTRALI1201002 1010107 13800) /01NP (3807119 
SMOZ © SmNsEhAd 
200 00 290 10] mar 
Khel 2O1 61 
e421 OK 
OlmP im) = ZAVECUCIKND/ROR, 
O1mPinmel) « SnoZ 
OlmPinmwe2) = MACH 
Keges 
29¢ CONTINUE 
OLMP Chee) aL AVES (WC CROMEGOT I /EXTRAL 1415 
300 CALL CCCIHCILDeLZEDD oMCHL DCC oRERD 
€ SET UP WOIZ ADM FOR AB 
IFIKERD 92003300310 
310 WRITE (60920) 
320 FORMAT (1H0 330X390 ERROR PROGRAM Co ALL CASES TERMINATED } 
CALL COREICE1) 68290CC) 
GO TO 90 
930 IF t oMGT. ARUM ) “0 10 360 
CSSSESS 6800000 0358066 8606808 FOOSSOEHDS SOL OFEEHSSE SOC OOS SEE OOEESOOVETOOS 
¢ MOVE OUTPUT FROM C INTO INPUT BLOCK TD & 
CECCETS F689 6 FS 0000000 6082 0000020000008 06 0490890 20AS 50866000 88640000008 


340 WP 8 CIS) 

twos 6 
00 350 T = le we 
two e twos 2 


ARITCS pe 221 1wo ’ 
AMIT(TO993 © Z2201wOo 1) 
AMIT(Te60r © ZZUtwOol0]) 
350 COMT TRUE 
€ ZEROES TO ERD TABLES 

AMIT(MPS 1) © 000 

AMTT(WPO3S}) © 620 

AMIT(MP61} © 0.0 

1f [| MACH olEs O00 ) U1Bam © T7203) 

GO 10 612 


360 IF ( .MOT, BRUM sAND. SMH EO. 0.0 ) 6 10 $20 
(0090000000000 0004056 050655 0ES0OFA5500 0686S SSSSS CESS OES SO EDOOLESCOOFOOES 


¢ MOVE OUTPUT FROM PROG C I'L ‘Wut SLOCK FOR PROG 
COC0002 0005009 0050005 505 SSEF 0000! 2000000 25ESSSE CSET ETS 


370 LO 380 1160200 


BOCCGES 0808 


9860 BitrezZzety 
Otape22Z1 2051 
€ 
309 CONTINUE 
HMRUM © (CH eGT. 000) oORe BRUM cANDe SUE GE ee} } cAMOo I CHsEQe000? 
e MOZILE ADMITTANCE TS INPUT TO PROBS As Be 


(OOO 000 5000006500020 PEHOHEROESESOSEORO SEC OCEE 
¢ 
400 1F ( .NOT, ARUN } 00 TO 080 
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#O6OCOSSEC ESOC OSES 


C1anl19@ 
CHAmi 200 
owniz36 
ont 220 
(A230 
CHa 200 
cnaniz3e 
Cran 268 
ca 270 
conn 269 
Ceam200 
cnani3e8 
crenni31@ 
um 326 
ea 393 
Omni pes 
Can} 390 


cham2110 
CHAMI120 
CHAN2130 
CHAM2Z140 
CHAMZ15U 
CHAM2160 
CHAMI170 
CHAMZ180 
CHANZ190 
CHAM2200 
CMAMZ210 
CHANI2720 
CHAMZ230 
CHAM2240 
CHAM2I290 
CHAMZ260 
CHAM2270 
CHAMN2260 
CHAMZ290 
Che M2900 
CHAM2$10 
CHAM2320 
CHAH2330 
CHAM2980 
CHAN23SO 
CHAM2360 
CHANZ970 
CHAM2380 
CHAN2590 


< CMAN2000 


e190 Cac LONGE Cmam2ei0 
€ CHam2&20 
1F (a €@ 1020200008706 CHAM2630 

O° MOTTE €608301 - CHAM2040 
= 77 FUGMATE LO eeOKeT Ae FSECR PROGRAY A) CoHAM2650 
TOLL COME 6 Batlle OCs Cad cuam2e60 

oe” COMTI ME CHAM24 70 
6% 10 110 cwanzeso 


COSC COOE OEE EOOOSOS OE SEETOSEEOEEOEESEEOETOLESEOTOSODCOPOSOSSESEOLOOCEDE CHAMIAYD 


. SET UP MAIN CONTRO, CATA FOR USE IN B Cran2300 
7 Od0 OOS OFS SSOCEESCCOEODED Seeesssesaaasenss CHAMm2§510 
e350 16 ¢ .mOT. Bum 2 GO TO 530 CmAM2520 
Btpseweeis CHAM23S 36 

St PecERTRAC Sr 7ERTRAL ILO) cuam2$eo0 
BIypeExTRariys - ChAN23380 
BtSpeEev@aties ChaM2560 

ar © OFaPey Fs CuHAN?S 70 

4 = DPMP E TAD cuam2500 

bad = INIT? cnan2s9o 
CIeCiny @6a) CMA" 2600 

DO 060 fel smOmEG CHANZ610 

PE Lote@roweite?s CHAmM2620 


e267 COnt int 
@rze7) © Oe. 41) 4 OP elas 
RED202 © DISter2s 


Atz7d © Ber CHAM2660 

DO 280 Ie 2. 20 CHAM2670 

TF IOUSTHLT 29090 069OeaTS CHAMZO80 

OPS BULOLOMMODISTLALIENT@aC ar CHANZ690 

At Ler SseDTSte dD cuan2700 

METI IOMETIICI,O cMaAN2 710 
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$3 © ROMEG-: 
e211 
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MWY te iese 
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CHL TOTO 6 BITTE LD. AmIT 2205 CMSCEKDs SLPRIRK! » 
CALL INFRA 6 BITES. OmITCe1 re OMCCER te ALFLICG? + 
reins 
2 28? fe le fe 
UF € YF 6Se St .OBe f oFSe 282 F 
eeeacTaet to} 
ALOR e TARE O71 
ALF) woFagt 1.31 
DENS (C1GA0-y} Bahoas ORE 2006721 BAR72OAL Pi moa Ple 
BOT © 1C 104806 1848-1 BARZ© ( A: ORL OA, ORE AL Pi meal Pls) 1/OERe 
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VOL TA52.00OMEGACES 
STewPelessaq routs 
COSPMTelOS 1 Owt) 
SlaveseSiic Tet TAs 
CTS TeeeCISe META 
CP 5 OORRE CCOSPm —BIMOST ura 
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ocho CjMAReA-Copje( ER 
CAP Jothl sGateareiuis Ca-C1 8a 
RORISTEGO”JoE OCR os C1 BAMOA-CODDO(EROLT o£ 1eDR) 
CONOR eFLOIMOCT-CeOts 
COM) oF MO BeCt-Ceott 
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CAP oF Te rger 
DE tmneca COPE OCAPL OC AM, 
CaPwel( Career alloc OCEP Zs CK ueN 
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MADE LS 2181992 To OMDEL 
G3 TC LPF 
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DEL YAO MMOFL OMEGA 
Wt -NIMENSTONSLIIFD MESATS 
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REtuaN 


End 


Page 29 


BtWGre Wie 
Pals FSTR Fre FL me ERIBSE IGT SOUND eC Fhe 


SB 2S> OCS °S FOO LOMGI TBI Me, MORE 1/2) STC ICRC) 


1M 060 
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LONG 990 
1ONG1000 
LONG1O10 
1ONGIO20 
LORG1030 
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CONS TOFOITELOL? 
xnm}Ne3300.0 

CALL PAGEC TO) 

CALL Ovche (ROOOFE) 
DO 26 1016133 

Atl) oF Ime) 
Btreo.0 
ERCTF IKI Omw) 

1h 0Onw 

TE CCF -99.O3 THe PH o 30 
CALL PAGE 14a} 
WRITE (6000; 

FORMAT 1340050 PROGRAM F 


ImPUT SOLVE FOR Mi) ANDO TAUCH) 


1/7 19% 0 OH OMEGA 1D 9X c OPT INT © LOX oGetesT IMT) 


00 $C I 91500 


WRITE (OeOOLOMGALTDOMTAIMTIS Ds MTTENTEND 


FORMATE 1d 610K 03F 1606) 
COMTI mue 
90 120 1 014% 


RMUSeLMTRIMTCL OMTRIMTIL) © MTIITMTCLICHTILMTAS DDSI 2. OOMTRIMTITDD 


OmOm o AMET) = MTRINTIND 


CALL QUAD 10MOM OHTTINT CR De TAULT DD 
TAUET DOL TAULLSCAC1 IORI. 399SS9I SEAL ZISOMGAT IN? 


yoreat2»eCO 

OME tT eOMGAIT) 

CALL OVCHY KOOOFR) 
GO TO 180290) eK COOFK 
MERSC 

GO TO 220 
TFCKMETI FLAN LAMINA 
TFCKMCT SEMIN) 110022001220 
Rrmelwennt ys 
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COmT taut 

DO 136 1e1+100 

FOuTtr) © BCT? 


IF ¢ CF-9.0 } 16001800180 


CALL PAGE (49) 
WRITE (60330) 


150 FORMATI1MO20H PROGRAM F OUTPUT //21XeTHECICPS 113K cOH1OMGAIDe 


160 
170 
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IRSReTHTAUIMS) o1GKo 1mm D 
DO 360 121240 
FCCPOSET OCOMSTOOMKIT) 
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FORMATI2IX F721 rLIK OF 1005 19K oF 1050501 OX OF 10051 


MTTIMT Catt © 9,4 
HTRIMT (ALT CO60 
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eYay 
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mTay 
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aTay 
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aTau 
nTau 
nt au 
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mTAU 
aTay 
nTAU 
nTau 
wtau 
nNTAU 
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wmvay 
aTay 
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Yau 
wTAU 
NTAU 
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CALL INTOLFCCOSEE DoMTRINTI LD OF CMI MGHTRE ) 
CALL FHT@IFCCPOSCLD OMTIEMTID) cE CMINGMTEL? 
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CALK QUAD DNOM HTT 1» TAUMING 
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WRITE (60210) 
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WRITE (hel POIHTRUIN: TAUMIN OMNI No FCMIN 
RETURN 

END 


Cc QUADR LIStemee 
SUBROUTINE QUAD (AeBe ANGLE) 
TF16) 10950+80 
TFCAY 29030000 
ROTATE « 361415927 
GO TO 110 
ANGLE © 457123890 
GO TO 120 
ROTATE © &o2831853 
GO TO 130 
TF (Ad 6007" 70 
ANGLE © 341419927 
GO TO 128A 
ANGLE © 020 
GO 70 120 
TFLAD 209990190 
ANGLE © 165707963 
GO TO 120 


100 ROTATE © 000 
110 ANGLE © ATANIB/A) © RO'ATE 
120 RETURN 
FAD 
SERFYC KORE LIST sme 
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RETURN 

END 
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SINCLUDE seCtoR 

SIGFTC Ing LIST mee Sms 
SURMOUTINE IM ICTR ms) 30 
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COMMON /PROLCGs LOGIKISB)e MEADIIZIe Silo SLPe CORs tas 30 
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EQUIVALENCE (LOGIRES2 cEM@ud « (LOGIKELO) sc Meme 2 o (LOGIKIPDIMUEIINS §©660 
LOGICAL LOGIae EMMe Bums Sile SLZs CORJe Rum ims 70 
DIMENSION 9 6EaH1 612) tay) = 80 
¢ LL  ) 
tf ¢ S19 GO v0 20 1%3 100 
90 10 f © te 0600 my 130 
a 30 EJDATACT? © 0.0 ims 120 
SL1 © eTRUE. tas 130 
$6 10 30 my 146 
¢ ms 130 
2° @EaD (139) Eymare ies 360 
GaccePace 33 ms 170 
70 CALL =AS298 { ESDATAs HEADI, ME 1 ims 1800 
° 1F ¢ ME OME 1D Caaa oexIT 1s 190 
WRITE (339) Cspata ms 195 
BACKSPACE 13 ms 396 
1F ¢ oWOT. JUN) GO 10 40 1" 200 
calu JJ ims 210 
40 1F ¢ Efum .OR. tMum d CALL INTIS ims 220 
€ my 230 
50 RETURN "ms 260 
eno In) 230 

sinere sect LIST omen SECT 
SUBROUTINE JIS JECT 10 
¢ RCT 20 
Cc eceseeseoese OfCK MODIFIED 20 AUG 67 EPPIEOE COD PODER STOOS JECT $0 
c ACT 40 
DIMENSION THETI(2090). R16169096 TeU6 10002 CT $0 
DIMUS 10m OFLE 10001 0X 11000) 261000) oNTVPEE! 1000) M0222 5ECT 60 
Leo TMEV Al 102d AS (20) oWTE (3000 1 AAMT E1008) CARP (1600) oMBAMDL LOGO eKMREI JECT 70 
21000) oXMUTOT! 20) oAXTV1200) AF TV 200) oXR6 1000) e VVC S000) JECT 80 
COmmgn /JECTORs DATA | 9600) SECT 90 
€ JECT 130 
EOULVALENCE (DATACSIORMD, (DATALOD oXND JECT 140 


EQUIVALENCE (DATAL3O22F oTHETID oc (OATAL2ZOTIDORIDCIOATACSOZSI oes § SECY 150 
EQUIVALENCE (NT VPEE cDATA( 200115 + CAXTY CDATAI10019) eC APTVeDATAIIZ01IJECT 160 


120 (COATACEO21) oRKd (DATA 20230 oVV De (DATALA921) oxRITOT) JECT 170 

¢ SECT 30C 
MERRORAD JECT 200 
SECT*OATALS) JECT 210 
WTODATALOSTS) JECT 220 
RIMJOOATALOST4) JECT 230 
XOMSPATALOSTS) JECT 240 
FECsDATAIOS 76) sECT 290 

> OF FCaDATAIOSTTD JECT 260 
MT SOAT AOS 70196.0001 JECT 270 

a POL aRedarales72) JFCT 280 
: NESDATA(9S7290.0001 JECT 290 
RONSOATALOSTO: JECT 300 
ROFeDATAI9SE0) weet 919 
Cmmeaxes.0 t 1 320 

« PRADJSE160 w#CT 3390 
PFADIE1 0 ECT 360 
COAeDAT AOS OS) JECY 390 
COPaDATAI OSG) JECT 960 
RFCuDATALOS9S) JECT $70 

Neer enaratesery JECT 980 

PFEC © DATAIOSNO1710000 SECT 990 

PREC © DATA198O01730000 JECT 400 

1F (KM FQ, Of 1 XM © PEN HeCT 810 

TF ( RM 6FOo Moo 1 UM © 1RD60 JECT 020 

+ weayo FCT 430 

CM 6 eo OV1OS37E0 1 Jt T @au 

v4 FRE9279920030090 1€C1 #90 

cre an do 1&1S026/5t6 7 JbCT OC 
oO 10 «0 JECT 70 

$s ROTATEQU,* JECT @80 

OO TEEPOL AR GTM o* HECT 490 
80 LEESEFT= 1 eDIAN eH ene JECT 300 
OC 10 7H Leteme Jecr $10 
Keone! JECT $20 

MELA THeDR TAIT R DO,0N0} JECT 330 
RIPPeMATSI Ke) JECT $40 
VEDPeDATAIR KO?) JECT 350 
NIVPEE LL ODATA(KK S314 20001 JECT $60 
SAVESSQRTERC TOKE TIOVETDOVUIDD JECT 370 
VUEPPOCATAMECVEED KOT) ) DOROTATE JECT 580 
RITVesAVE JECT $90 
Kees CT HOC 

: 70 CONTINUE JFCT 610 
GO 70 270 JFCT 620 

a4 00 170 Tsdenr JECT 630 
Kets) FCT 640 

MELE CT eDATAIKK 1*.000) JECT 699 

NCP sDATACK KOT) JECT 660 

“2 r YUL sDATA(KKe2) JECT 670 
NTVPFE CT ISDATA(KK 031600002 ICT 680 

¢ THE FOLLOWING SEGMENT OF CODING REPLACES SUBROUTINE MOVKYoes JECT 680 
weoeevety op 90+1900140 JECT 700 

90 eet mete 1009310s120 JECT 710 
tan ROTATE © 441615927 OT TA 
GO 10 160 ECT 790 

die ETE © 260f=25 JECT 740 
120 ROTATE © 6,2831853 JECT 7490 
GO TO 164 FCT AN 

1A ROTATE © 050 JfcT 770 
GO TO 1860 JECT 780 

140 te OC utte yd 1000150130 JECT 790 
150 KET) © § 608-35 1€CT fO0 
ROTATE » O00 JtCT 810 

160 SAVE @ SORT C REDDORCTDOvEtpevetE 3b JECT 620 
SAVECSORTCMELIOMCTDOV¢LIeVENDD JECT 830 

YELP STATAMEVEEIZRED DD PeROTATE JECT Ban 

KUT beSAvE yFey as 
Ceges Wet owia 

176 CONT TNUF CY BIC 
GO TO 274 wet ABO 
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189 1FSSECT<1569190 03900200 JECT 890 
190 MOTATE*O.D SECT 9090 
200 00 219 1st omE JECT 910 
Rieke} JECT 920 
MELE CT PeDATAIRK 6.0001 JECT 930 
REPDOPATALEKOS) SECT 900 
VOTE eIDATAIREO210.01 70992 1OROTATE JECT 950 
“TVPEE IL ICOATAIKK 23 14.0001 JECT 960 
cones JECT 979 

219 CONT SnuE JECT 980 
CCRC CCCC ECR C CEE CE CCC EC CEC CECE CCL CECE ECEC CECE CEE CCC CECLECCCCOC CCE CC CECE SECT 990 
c CALCULATING AREAS OF ELEMENTS JECTION 
P CCC EEE CEC EL CEC EECCCECE CECE ECC ECE EEC ECE CCC CECE CECE CEE CECE CECE CEL EC CECE ECT 0° 0 
223 aMnetee SECTILNO 
30 260 fat ont JECTIOSO 
KRowke} JECTIO“O 
NSTeDATAIKK 4.0003 JECTION 

RK PORK SECTIO“O 
JJonnl JECTIOTO 
NESDATAIKK 1120003 SECTIOCO 
ALTYINST 10000 SECTIONS 

IF 197952900730 JECTII00 

23° DO 260 Jeol oNF JECTINIO 
SJORK LOI SECTIIZG 
ARTYVENST DOARTVINST 101 o TOSSO9I OUDATAL ISI ODATALSIIND JECTIIIO 

240 CONTINUE JECTIISO 
250 NMleJ Je. JECTIISO 
tos 300s), JSECTIIGO 

NF eDATACINII+.0003 SECTIITO 
AFTVINST 100.0 JECTIIGO 

TF CMF 126002800260 JECTIISO 

260 390 270 JFel NF JECTIZ00 
NNeNN) + IF JECTIZIO 
METYENST DeAFTYINST 166 5 78999810 ( DATA CMD CDATA LIM) DD JSECTI220 

279 CONT SUE JECTIZ90 
280 CONTINUE JFCTI2Z40 
ANTOT#O.0 JECTIZSO 
AFTOT#O.0 JECTI2EO 

00 26 101 ME JECTIZ7O 
NMOMTYPEE CID JECTIZE0 
AAXCLVOCAXTY (0) JECTIZ90 

ARF TVOAFTYCNMMY JECTIION 
AXTOTSARTOTOAAK(T) JECTIIIO 

AF YOT=AFTOTCAAF II) JECTII20 

290 CONTI NUE JECTISIO 
ANTOTSARTOTSSECT JECTIIOO 
AFTOTOAFTOTOSECT JECTIISO 

bad a ddd ddddddadddddddddddddddaddddddddddddddddddddducddddeddddddudadadae ast) |) 
¢€ CALCULATING RADII AND ANGLE GOUNDRIES JECTISTO 
POCOCOCELE COE CO EO COE LOC CCE OOO CLEEE COE LOO FE CECOO EE COE ECOC EC COE CE EC OOCOCE JECT1 980 
300 NPTSe100.0/SECT JECTI39O 
TPSeNPTS JECTICOO 

DEL THO THENL /TPS JECTIAIO 
AREA © 3.141592 © Rinjoe? JECTIAZO 
ASECT © AREA/( KMOXNY JECTIAIO 
R11) © 000 JSECTIO’O 
ANP © ASECT © 18000 / 36141592 JECTI“GO 

DO 310 Tf © 2021 JECTIOGO 
RIT) © SORT ( ANP © AI T-31002) JECTI“TO 
310 CONTINUE JECTIAOO 
R122) © 0.0 SEC TIAGO 
XnuMe0 20 JECTI500 
NPTS1 CPTSe) JECTISIO 

00 320 Ja2enPTS) JECTIS20 
RUM 0 KN 160 JECTISIO 
THETAL JI ODELTHOKNUM JECTISAO 

320 CONTINUE JECTISSO 
SCCCCCOCECECCCERE CCE CCE LECCE CECE ECCCC CECE CEC CECE CCCC CECE CEECCCECCCC LECCE JECT ISGO 
¢ CALCULATING WEIGHT FLOW JECTISTO 
COCECCECCCECCCCEE CEC CECE CE CEE CE COCO CCR CECC CEE CEC CE COE ECCCCCCCCCCCOCCCECC JECT 3980 
WFTLOWT/( XMRO260) JECTISGO 

WET L ONT=WFTT JFCT1“OO 
wET20WFT] JECTIOIO 
WXTQEWKT] JECTIO20 

TF (OF FC 134003600330 JSECTIOIO 

330 WET 2EWF TLE eO-PFFC? JSECTIO“O 
360 IF UPXFC 136003600380 JECTIOSO 
350 WXT2eUKT1E( 1 60-PXFC) JECTIOGO 
360 SIGMENCWTs 1 40161992ERINJORINI) JECTIO7O 
370 AXTMAX@AXTOY JECTIO“O 
AXTMINSAKTOT JECTIO9O 

AF YMAXSAFTOY JECTI700 

AF IMINSAFTOT JFCTITIO 

380 AFFC#0.0 JECTIT20 
ANFCen.0 JECTITIO 
TFLDFFC 140028005390 JECTITGO 

390 AFFCAFECO, TESSOReDFFCODFFC JECTI7SO 
AFTMINSAFTMINOAFFC JECTI76O 

AF YMAX CAF TMAX~AFFC JECTII7O 

400 IFLOXFC} 42006200410 JSECTI780 
910 AXFCUXKFCO@, 185398 GDKFCODKFC JSECTII9O 
AXTHIN@AXTHMINOAKFC JECTIO00 
ANTMAR GAN TMAX©ARFC JECTISIO 
420 AXMOMs (AXTMAKOAXTMIND /200 JECTIOZ0 
AFNOMs ( AFTMANGAFTMINI /200 JECT1O30 

TF (DFFCODKFC 1400 08600490 JETTIA‘O 

430 WXT26 (WHT 25 AKNOM) /¢AXHOMOAXFC) JECTIBSO 
WET 201 WFT2 AP NOMI /(AFNOMOAFFC) JECTIOGO 
SON XMR] OWXT2/WFT2 JECTIOTO 
490 ETOF 80.0 JEcTI8aO 
ETFF e060 JECTIO9O 
QXewxT2/AaxTOT JECTIVOO 
QFewFT2/aF TOY JECTIV’O 
CALL OVCHK {K000FX) JECTIVNZO 
460 00 660 191 NF JECTIOGO 
NBANDI 100 JECTIOSO 
WFELeaar (1) eOF JSECTIOGO 
WHE] eAAK( 1} 0OX JECTIOTO 
WELT DewFE LOWKEL JECTIOND 
XMRE TT bewxE 1 /WFE) JEC T1990 
ETOFSETOF oWKE1 JECT2000 
ETFROEYPFOWFE) JECT2010 
CALL OVCHK (KOOOFX) JECT2020 

GO TO (4702480) KOOOFX JECT20I9 

470 NBAND( Tbe-) JECT20e0 
480 CONTINUE JECT2OIO 
ETOF eSECTOETOF JECT2060 

FTF OSECTOETEF JECT2070 
WHREL MOE TOF /ETFF JECT2080 
AXXTTCARTOTOARFC JECT20%0 
AFFTTSARTOTOAFFC JECT2100 


Mt ES Se arin en re ketene Mee a 


wet2 © Eee SECT2210 
wut? © E1Gr MCIZI20 
WT owFT2/ 10 0-PF FCS ACTZIU3IO 
WRT 1 eWI2/1300-P2FCD JECTZ1460 
ETS © eT OAFFTT Ar tOT FECT2I90 
WHT} © WRILOARETT/ARTOT JECI2Z160 
SMAI A JOwNT 1 OFT. SECTIUTO 
Deo C(UNTE/ICOROARRTY 25 007302,2960266 JECT2100 
wtOTeunTi courts RCTIIVO 
OPxieD/80K SCI2200 
‘* De CIWFTLsECOFORERTY 1 9002107.2360208 HCI2Z10 
DPF 1 =D /ROF JEC12220 
TOFCF © wETBOII.0 ~ OXECS JECTIZ90 
TOPCH © CTOPCROIARETT - ARFC}I/ANRTT JECTZ240 
TORCE © WET? = TOFCE HECTZZ50 
YEFC = WT}01200 - PFFCI HECT2260 
IEEE © IYEFCOLAFETT = APEC D/AFETT SCTZZT0 
TREC «© wT) =< TEPC HCT2260 
° VIM JLCSORT(192732600PK) /ROKI) HCT2290 
VIMIFCSQRTI69273.600PF1 /ROF 1) JECT2300 
RRvatP: IGMEMOSECT JECT2310 
O90 WRITE (608201 HACTIZ320 
GO TO $00 HCIZ330 
$00 WRITE (65760) ANTOT CAF TOT AKFC oAFFC GALT! cAFFIT DPX) DPF} JECT239400 
URITE CheTTOPUTOT cPMELMeXMRIMIoFTOR cETEF STORCE STFEC. UET1e WET] oJECT2950 
IMIR oVIN IE JECTZ360 
PF FC ePFFCO1O060 SECIZ370 
PREC OPXFCO 10000 SECT2300 
$10 WRITE SO sal OMT eE SCOR JECTI990 
WRITE COeTROICOF oPEFC oPKEC oOFFC COXFC oF FC oRFC ROX OROF JECTZ400 
CALL PAGE 70) JECTZ410 
WRITE 160810) JECT2020 
WPGE C49 SECT2OIO 
DD $40 Tel NE JECTZO“N 
IF CRPGE~2 392005200530 JEST2ASO 
$20 CAL. PAGEt TO FECTZ460 
NPGE SHPGE O08 JECT2470 
JECT20BO 
$30 JECT2490 
Bextlsesimevisys JECT23900 
MeKCTICCOSEVEITDD JECTZIS10 
WRITE (GcHOODMELE( TP oNTYPEE( 1b oX (1) DEGREE Ac JECTZ2920 
340 CONTINUE JECTZ930 
CALL PAGEL TO) JECT 2540 
NPGE 946 JECTZ950 
JCONT SO JECT2960 
WRITE (6.000) JECTZSTO 
WNCATES JECTZ990 
DO 670 featont JECT2990 
RKONNS] JECT2600 
NSTEDATACKK 140000) 3ECT2610 
KK PSKK OL JECT2620 
JJORR) JECT2630 
MXCDATAIKCKI 1400001 JECT2040 
JSCONT SJCONT+NXS2 JECT26SO 
IF (NPGE- CONT 155005500560 JECT2660 
$50 CALL PAGEL 70) JELT2670 
. WRITE (60860) JECT2660 
JCONT SO JECTZE90 
560 AXMAK®000 JSECT2700 
AFMAK#000 SECTZ710 
WFE1#000 JECTZ720 
WHEL Ce 4ECT2730 
WRITE (60BS0INST OMX JECTZT4O 
UF 64K 197005900370 JECTZ750 
7 570 DO 580 JXe1 NK JECTZ760 
JIOKK IF IX JECTZ770 
Ate TESIOPIEDATALJJIODATAC JIT JECTZ760 
Beaegn JECT2790 
WRITE (He RGOIDATAL JID OAC JSECTZ2800 
AXMAXCAKMANGA JECTZE10 
WXEL OWREL OR JECTZE20 
580 CONTINUE JECT28390 
$90 KNLeJJ+ JECT 2000 
MNORNT JECT250 
NF eDATA(NNY 3400001 JECT2060 
JCONT # JCONTONF 42 JECTZ870 
TF (NPGE-JCONT 1600106000610 JECT2E80 
600 CALL PAGE (70! JECTZ690 
WRITE (60800) JECT2900 
JSCONT ED JECTZ910 
610 WRITE (608 TOINE JECTZ920 
TFCNF 362006400620 JECTZ9390 
620 00 630 Jfs) NF JECT20€N 
MN ONN] OIF JECTZ950 
Ame TOSIODLOOATAINNIODATAINND JECT2960 
Beacegr JECTZ9TO 
WRITE (6 eBSOIDATAINND cAnB JECT2¢000 
APMAXSAFMAXOA JECTZ99O 
WFELAWFE]OR JECTI000 
630 CONTINUE JECTIOIO 
$40 WIOTLOWXELSWFE) JECTIO2O 
TF (WEE) 166906500669 JFCTIO3O 
O50 WRITE (60900) AXMAXeAFMAR pWwIOT} JECTIOSO 
GO 10 670 JECTIOSO 
660 RMR OWXEL /WEF I JFCTIOGO 
si WRITE (60890 I ARMAY AFMAX sWYOT] oKMR} JECTIO’O 
670 CONTINUE JECTIOBO 
680 READ (13) DATA JECTIOVO 
BACKSPACE 13 JEC T3100 
LINKNT © 0 JECTIIIO 
ASS © XMOLN/WT JECTSIZ0 
DO 700 J © 1seNE JECTIIIO 
é TuUtJ) © WTEC JIPAIS JECTILGO 
TF CLINENT eGTe 63 GO TO $90 JECTIISO 
CALL PAGE 70) JECTIIOO 
WRITS (69910) JECTIIIO 
WRITE (65930), JFCTIIBO 
LINKNT © $0 JFCTILEO 
$90 WRITE160920) Je Kidde Veshe TMUL I, JECTI200 
RUMENT © LINKNT © 2 JECTI210 
700 CONTINUE JECTIZ20 
DO 705 J = le NE SECTIZ22 
RRCfe os MEED JFCTI226 
Tos yytt) © ved) JECTS225 
WRITE (19) DATA JECTIZ30 
AACKSPACE 13 JFCTAPEO 
wMyeared.f JFCTAISO 
AIS © AISOSECT/ EW IFCTI26O 
BDO 4A 301020 IFETIZID 
TOTWaN, A HEC T3240 
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REX © ERY © ASECT SECTIZ00 


Jes 3ECT33990 
00 710 Let omt JECT33910 
TF CREED chE. RI JOD) cAMbe KES) oGTe MIJS2) TOT » TOTW © wTECEs §«©JECT3S20 
120 CONTINUE SECTIIIO 
RmyTOTIJ) © FOTWOATS JECTI940 
IF ¢RUUTOT ¢ 5} -KOLMAK 173007300720 JECTIISO 
720 XeMaxexmyTOTs J? JECTS960 
730 CONTINUE SECTSITO 
CALL PAGES 70) SECTS3B0 
WRITE 167901 JECTS39O 
00 740 Jote20 JECTI£OO 
ARAMARK OXMUTOT I J1/ UAE JECTIONO 
WRITE (hoT- CIR Jol) oKMUTOTI I) oXXMAK JECTH420 
7a0 CONTINUE JECTI4I0 
RETURN SECTIONS 
TSO FOMMATIIZE to DolK of To30} 3X eF Gee) JECTIOSO 


TOO FORMATI/ 4. £028rthesccPROPELLENT ORF ICE AREAS C// cOXeSOMELEMENT TOTASECTIOOO 
WL ONIOIZER shea oF 2168 0BH SQ. Ie olSKeSOMELEMENT TOTAL FUEL JECT3470 
2AREA 2sFL 2B eB SOo IMoo/ ONO IOHTOTAL OXIDIZER FILM COOLING AJECT3O80 
MEA CoF Le 2 eBH S2o IVeeISKSSOMTOTAL FUEL FILY CCOLING AREA oof 218 JECTIOVO 
SoSH Sle Theos oP oIOHINJECTOR TOTAL CRIDIZER AREA oF 310858 SQOJECTISOO 
Se Tool SNe IOHINITCTOR TOTAL FLFL AREA SeF 11 BsO4 Sy IMo ess oSKIECTISIO 
CoSEMBoveeINJECTOR PRESSURE DROPS FOR ABOVE INJECTOR DESIGNe//e9Ks JECTIS20 
TZAMORIDIZEP PRESSURE DROP wofbod eIHPEI oISKeZOHFLEL PRESSURE DROP © JECTISIO 
BoP bel obht PSbes// SNe TIM ccecPROPELLANT FLOWS AND INJECTOR VELOCIT I JECT3540 
ES FOR ECVE INJECTO® OESIGNOs? JECT3590 

TT0 FORMATIZ oh OXeLOHTOTAL WEIGHT FLOW cof bolo TH LR/SECo/ 042 oSOHAVERAGECTISOO 
JE MIXTURE RATIO OF THE ELEMENTS oof To des eG 1ReSOHOVERALL MIXTURE RAJECTIS7O 
2TIO FOR THE INJECTOR 9 oF 7030// 09K e SAMELEMENT TOTAL OXIDIRER FLOW JECT3S8O 
3 BeFGeleTh LO/SEC eZ. Xe2GHELEMENT TOTAL FUEL FLOW 2 eFbelo TH LB/JECTIS9O 
SSEC oso IN SSMTOTAL OMIDIZER FILM COOLING FLOW = eFboleTH LB/SEC e23X oe JECTION 
S2OHFOVAL FUEL FILM COOLING soF ele TMH LB/SEC 7 s9X o3QHEMZECTOR TOTAJECTIONO 
OL OXIDIZER Ficw oe FbeleTh LO/SECe2IXe2SHINJECTOR TOTAL FUEL FLIECTSOZ0 
TOW oF Tole TH LAEEC 07. PKs 3STHOFIDIZER OVERALL INJECTION VELOCITY © oF JECTIO30 
belo Th FT/SEC 1OXo33MFUEL OVERALL INJECTION VELOCITY Sof bole lH FT /JECTIOSO 
OSECos/5 JECTIO‘O 

780 FORMATIOXsZ3HFUEL LOSS COEFFICIENT =oF S030/e9X s2THPERCENT FUEL FIL JECTIOO0 
IM COOLING =sFSele/seIReDIMPERCENT OXIDIZER FILM COOLING © oF Sol0/ 9X JECTIO7O 
ZeIOMOSAMETER OF FUEL FIL COOLING ORFICE seF TeSe7SH IMe (NOTE eo THE JECTIOSO 
3$ MIGHT BE AN EQUIVALENT GIAMETER FOR MULTIPLE-ROw COOLING) + /09Xo JECTIE9O 
SLZHOLAMETER OF OXIOIZER FILM COOLING ORFICE oof ToSeZ3M [Me (SEE ASIECTITON 
DOVE MOTE) o/ 69XeSIMNUMBER OF FUEL FILM COOLING ORIFICES PER INJECTOJECTSTIO 
GR 9 oF 503 0/eIXeS4HNUMBER OF OXIOIZER FILM COOLING ORFICES PER INJEC JECTSI2Z0 
TTOR soFSe0e/eIRelGHOKIDIZER DENSITY oF TeZ ehh PCF s/s9XeTOHFUEL DENJECTS730 
BSITY seFGoleel PCF.s7/) FECTITSO 

190 FORMATI//// 03K oS 1GHSECTION SeeRESULTANT FLOW DISTRIBUTION (My) AS JECTSTSO 
1A FUNCTION OF ONLY THE RADIUSs ‘Eee AVERAGE OF MU IM A RADTAL BANJECT3760 
20s PIF eSZAKoGHRADIUS 15K o SHIR 11K ok SHMUERD /MUIMAXD 9 / o FOX oSMIMJECTIIIO 
deols) JECTSTOO 

BNO FORMAT(13X0139016X012 4X oh1 10KeF 10050) SECTIT90 

G10 FORMATI// 22K. GONSECTION Zecee ELEMENT LOCATION AND INJECTION TYPEJECT3B&00 
Le/SL/ oY Xe THELEMENT o15Ne GHTYPE SL IK olHRol Ie SHTMET Ae 1 IN olHKoIOK eI JECTISIO 
2MYo/e13Xe BHMOoeLEKe THMOeesOXe SHIINCHES) 012Ks BHIDEGREEDeI2Ke OHJECTIOZO 


DUINCHES) o 12K eBMIINCHES} 0/7) JECTIBIO 
020 FOAMATIZ// 93K TOHSECTION Leos oMISCELLANEOUS INFORMATION FOR INJEC JECTIBRO 
1TOR DESIGNED BY PROJECTS af? JECTIOSO 


030 FORMATI//e5Xe &6MDeccoec INPUT INFORMATION USED IN COMPUTATIONS o/JECTISOO 
T/0Me JIHTOTAL CROPELLENT FLOW eof bole TH LE/SEC0/ 09K SPHTOTAL NUMJECTIOTO 
2BER OF ELEMENT TYPES (SYMMETRICAL SECTION ONLY) © 0330/99Xs SOHTOTIECTS: 
SAL NUMBER OF ELEYENTS (SYMMETRICAL SECTION ONLY) © si do/09Xe 27HOKJECTIEIO 
SIOSZER LOSS COEFFICIENT sFSe30) JECTS9NO 

BOO FORMATI//02Ke TOHSECTION JoooeTYPE DESCRIPTIONSORIFICE AREAs PROPEJECTS910 
ALLANT FLOWs AND MIXTUSE PATION// 03X09 QAHTYPE — Omman—omm OXIDIZER JECTIOZO 
2ORIFICE DATA ----~00------- o> FUEL ORIFICE DATA m-=--2— 00/7 o93XeJECTIOIO 
3 SHTOTALeSKe SHTOTAL ©5Xe SHTOTAL 9GX0 THMIXTURE 0/0.2Ke GHNUMBERs 34XJECTIOS0 
Go GHNUMBER 933K BHOXIDIZERSSXe SMFUEL 63X01 MPROPELLENT 2 3Ke SHRATIQJECTS95O 
Se/olaXe2HOF oSXe SHOIAMETERAXs GHAREASGXs SHFLOWs 7X e2HOF 05K» BHDIAJECTIO“O 
OMETER 4K SHAREAGE Le SHFLOWSOKs SHAREACONs GHAREASGXs SHFLOW RATES JECTIO7O 
T/eLAKe THORIFICEs6xs3HINeeSXe THSOe INeeeXe GHLB/SECO3X» BHORIFICEJECTIOSO 
BS 95K SHINS aXe THSDe INDOAKs GHLB/EECOIKs THSQe INoe3Xe THSQe INeJECTS990 


Geeks GHLB/SECe// 0) JECTROOO 
850 FORMAT(3X91357Kel2) JECTAOIO 
060 FORMATI22% oF 6sGe 3K oF Bebe 2K oF Bob) SEC TAO20 
870 FORMAT Sex5121 JECTAOSO 


6) JECTAOSO 
&eIKeF Toad JECTAOSO 


080 FORMAT i62x~ 
890 FORMATI9IX,+ 


FOO FORMAT LOIN F605 02X oF SoS e2X oF Oe Se2KeOHINFINITY) JtCT e060 
910 FORMAT (//‘JOXISHELEMENT RESULTS 1 JECTAOTO 
920 FORMAT (3PKTSOLINE DOs Se LORFUG Mol IXFINS 1 JECTAOEO 
930 FORMAT (//31XTHELEMENT1 SXOHRADIUS ISXSHAMGLERIXIZHDISTRIBUTION/9OX JECTAO9O 
LR LHCOEFFICTENT/36X3HNOe ISK THRADIANS 15X2HMU// } JSECTAIOO 
END JECTHIIO 
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Sterte Inses | LIST{NOS TY.) 
SUOROUTINE Insts insp 10 

c 1uso 20 
C Seececcceooso OfCK MOBIFIED 20 aus 47 SOTOSOSSETOOSSECCCCTE {Rm 30 
rd msn 40 
REM 19% 010K PT INTEL IPZ tinge 30 
INTEGER OCHO. TIME 1m 60 
LOGICAL LOSIKs Mle Side EORJ» 1mm tas 70 
CONOR /PRELOG/ LOGIKISO)> SLie Site CORY 1nd 80 
Common s3ECTORs DATA ( 96001 IND 90 
COOMION =» GANA, Me WL e AVMs «OvMe Cum, Cull, Ch, Ct taun 110 
Blmgns 160 14s 120 

1 8 (100020 AVE 6 ¢ (3Ole GU Sere CVMR ¢ 3020 8 «18ND 190 

2 «CVE F900 MR «= C100RDs «THATA (3000). FIRST ¢ 3000 IND 140 

3 SECOmmt 3620 2 (2% wt (88te 18D 180 

a 19715920 CPPLSODe LPRISE)s CPMISOLe IPTISOD» OPTISOD: 1a 160 

s 110089. ¥ 130009 Igo 370 

7 CQUIVALENCE ¢ Bums LOGIKI9) 16 Ing 100 
1 CES eOATALL ID (um o@ATAIDID 18JB 190 

2 (mh DATALOD tv sOATAIODDoCSYN «= sDATAIBID» INJD 200 

S$ (POC oDATALI99D» 623 eOATACISDD» INJO 2130 

s (MMe de DATALIOZII Dol TMATAs Vo OATALZOZ1LII® 1N® 220 

S tu eBATACIO2ID Ds 18 230 

. 6 C10P =ePATAIZO} de imo 240 
7 (OPP DATALTODDcE2OR eDATACIZ0D 36 (OPM sDATAII700)6 1mJn 240 

e CIPT eDATAL22009c1OPT eDATALZTOD Del VFLPDATAI9S06) 16 1" 260 

9 ETFLR DATAL OSS?) 0 oC TELT COATAL9SOG29 Heed 

€ LP 0 


10 FORMAT (//12XDOURESULTS OF DESCRIGING FUNCTION // Z29K1OMELEMERT  INID 290 
23 GOMADIUSOKSMANELESZIONPRACTIONAL /19XSHONEGATXINO. TEIN =«=6—- TASWRADINID 300 
2 SROFLOU-RATELONZHFPLANINFRIZKINET 4/7 9 tnyo $10 
20 FORMAT LIOKF9 6 o1SoF 12 oF oF IOGoF IIS OBNIF IOS F 1NID 320 
30 FORMAT (77 12RE2HRESULTS OF INJECTION BISTAIOUTION EFFECTS 17 IND 330 
LOPSRIMOMEGCASK SMAVNGESCBUNTZSNCVNTINSHCV / SOXSIGHOMMEAL 1 o6KGHMIMAG/ NID 300 


27) 1njy 990 
40 FORMAT [99N5F10.4 3 Tayo 360 
S@ FORMAT ( GOXSHALL 471004 1 {nso 970 


60 FORMAT ( // SKSOHLMPUT TO IRJECTION OLSTRIGUTION PROGRAM// 1ONOHCONINID 360 
SSTANTS//1ORIOMUMBER OF OMEGAS © 139//1OK2Q0NMEDER OF ELEMENTS 215,1NID 39, 
3 190 FOR EACH OF See 20H SYMMETRIC SECTIONS, 1M. 392 
2/AUAAIZMRADIAL OLVISIONSIAM) © F506// 1ONZEMANGULAR OLVISIONS (XH) [NID 400 
9FSeO//1OKZTMACOUSTIC MODE MUMBERISYN) OF 7.4//24XSOHOROER OF BESSEINID S10 
GL PURCTIONSIV) © FS.O//1OK1THINJECTER RADIUS OFO.SeSiHe Ios IORSZHINID 620 
SRATIO OF SPECIFIC WEATSI(GAMMA) OF 7.4/ /1QNSTOMAXIINM PRESSURE AMPLIINID 430 
OTUDE RATIO‘POO) =F7e9///14KSONTRANSPER PUNCTIONS FOR LINEAR OPERATINIO 407 
TION//2OR1LGHPRESSUREC TELE) OF 7o3//2ONZSURADIAL VELOCITYITFLA) oF7.3INID O50 


O//ZONZTHTANGERTIAL VELOCITVITFLT) oF 709 1 1NID 066 
70 PORMAT (7/ LOXITHINPUT FREQUENCIES ) 1HJO sto 
80 FORMAT (/ 7X SF28,4 ) 1NJO 480 
90 FORMAT (///10NOMELENENT INFORMATION } TMD 490 
100 FORMAT (//DIRTHELEMENT LOXGHMADIUS ISXSMANGLELIZIZNOISTRIGUTICONs90OX I1NJO $00 

VALMCORFE 1C LENT / FOND. 1 TAIT SSRPMRADIANSISRZMU// 4 IMJ Sto 
120 PORMAT 622X550 132F 10090 1OXF 16.401 3NF 1006 7 Ingo $20 


120 FORMAT (1H3//26H TABULAR MOMLINEAR EFFECTS /11XQHPRESSURELIZ1 OHCOMTNID $30 
IOUST SON} 2KEHMADIAL 12K 1 OHCOMBUS TION 1 OX LONTANSERT IAL IORI ONCOMBUSTIONIN 1 S40 
Z/DIROMVELOCITVSZROHVELOCITYs OX SCL2MPERTURBATIONISAGHGAINIZX) > IhIN 550 


390 FORMAT ( GIF IG.902") 1 THIN SA0 
140 FORMAT (//99RESHOCCERCECCCS THESE VALUES PERTAIN TO A STANDING MOTRIN 57° 7 
, 1fZ eesecevese » IMsD SAS 
190 FORMAT (//S9KESHCCeReReCCse THESE VALUES PERTAIN TO A SPINNING ENIN ID $90 
10€ seeececess | 1mJ0 600 
3160 FORMAT (29KS2HTHE VALUE OF I2Z1T 1S ZERO AND MOY ALLOWED! thn 610 
P 170 FORMAT 2SKOSHVALUE OF BESSEL ARGUMENT TOO HIGH OR LOW Z = FIte&) INN 620 
1NJD 630 
100 OSCREe © 3 TN 30 640 
. 1F 6 MOT. TRUM} @ TO 390 min res 
CE © 100.0 THIET 66. 
oscrn © 2? inn 679 
60 TO 206 "NID 680 
100 Wwe 3 tn 490 
WEIL) © 160 Im tO 
200 Ting © 1 IN IM TIO 
Rims = OATAIOST4> tN IM 720 
WE © DATAI9S71) 2020003 INJN 730 
WS © DATA(S? + 02000) TNJD 73) 
WUNBR © DATAI9S99) ¢ 0.0001 TNO 740 
{22t7 © pavat2y tho 780 
K © ¥ 40.0003 tho 760 
KK eK oe} INJD 779 . 
b= 7eR- 1 TRIO 780 
Mele? tan Ten 
None? THIN 810 
IF ¢ ER ehEe 060 1 €1 * 60-0001 NIN 0 
tf ¢ MUR oLEe 0? nuMeR © 10 NIN PY 
IF ( CE oLEs 99.0 1 GO TO 260 (asf 650 
CALL PAGE { 701 1NJD 640 
WRETE (6060) MeteME WS oXMoXNoSVNe VOR INS GAMMA cPOOSTFLROTFLRGTFLT INJD 880 
GO TO 2200210)» DOSCRB TNO 860 
210 WRITE (6070) IwID 070 
WRITE (6080) (WCC tW)e Iwo Lene) tNID 860 
220 Linear e 6 tHiM 890 
BO 240 J © 1NE tNJD 900 
TP {LTMKNT 667. 09 GO TO 230 INJ0 910 
Py Catt PAGE ¢ 70) THI 920 
WRITE (6090) TN In O16 
WRITE (60390) IN 1D 960 
LiMn? © 90 tNJD 990 
290 WRITE (601109 Jo MRIS), THATAL JI. UIs) In 980 
LIMmNY © LIRKNY = 1 TNu0 970 
200 CONTINUE TNJG 980 
> GO TO 1260025015 OSCRE TN In 990 
290 WRITE (60120) 16 101000 


WRETE 160290) TIPPETT OPRI To LPRITD OPRETD STOTT OPTI LD ofedeShy = TNUNIIIG 
(C0 9900 90 0 000 90 OO OE TN I01020 
c RUNNING DESCRIBING FUNCTION IMPLIES THAT THE EXPANS!ON TH JOLAIO 
C — CORFPICIENTS Avie GVM» AND CVN ARE FUNCTIONS OF THE FREQUENCY SOMEGASTNIDIO6N 
¢ THEREFORE, AN OMEGA LOGP MUST BE ESTABLISHED, OTHERWISE, THIS LOOP IN ININe 


C 1S GONE THROUGH ONLY ONCE. IWnines 
CRAMMAAAMARARAAA ANAK AAA RRA AAA AA RAK CA AAR RAR AKRAM KER RRERE VOL ATO 
200 00 870 Iw se 1s Mw Inyni das 
WER © wCtiw) tw N1Ng0 

VOO = POO/ (GAMMAeWCA) tTaopy100 

end e rie Th 303110 

© 0. WNT 120 

BURT © O00 INJO1130 
LImEAT © © ImIDI160 

00 490 J = ts Ne TNIDEYSO 

Roe MAR(IIMING TS 101160 

OVA © ULsIsiReonny TH IHLi 70 

€ 0OG8eSCe8: SSOHODEROOSHEASESLOESSOOTEESHLESTECOSOSEESSSUSA EZ ECEDLEROTS TN N1180 
¢ TWE CALL YO @ESSEL REQUIRES THE DEFINITION OF THE REAL PORTION TN IND 390 
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OF THE ARGUMENT. ZI}te AMO THE IMAGINARY PORTION: 2126 FOR A 31%3D1200 
GIVEN OROER Vo THE OUTPUT OF BESSEL GIVES JOMIZ15¢ JOL(Z2)0 JIRIZIe IMIOIZIO 
JUNZ200 cece SUMIZID® SWItLZ}e THEREFORE s THE ANSWERS WILL BE 18301220 
STORED 1% THE FOLLOWING FASHIONe JiV-1) © FIRSTELis Jt¥) © 1m J01230 
FURSTINGe AND JIVOl) © FIRSTINI. 1%902240 
090000 00000000000 090000000009 000005 5050050 9800005006000 6 0000008 4600081805250 
wean e 0 1"JD1260 
Zt) © SumER SN ID1270 
2171 © Of 14301260 
CALL BESSEL € FIRSTELDs SECONDE1De KKo Zilts KERR } 1.1290 

1F 4 weme 5 27002700930 1”JD1300 
270 Sivu © FiasTine 19301310 
OSivm © f VEFIASTIoMn = ZILIOFERSTING 3/R $0JO3320 
Zeid © Sue 1401330 
212) © Oof 1451340 
CALL BESSEL f FIMSTIL) > SECONDIT)> KKe Eilte KERR 1 39JB13980 

Iv ( KERR) 260 0286.990 14JB2360 
260 GO TO (290032610 TINE 1*J5O8370 
290 1F CK Dt 300 03000310 $43013980 
300 © © & FIASTIMIOFIASTIND  FIRSTILIOPINETILD 2 / 200 $msD1300 
60 VO 320 $"901400 
310 D © C PIRSTIMIOFIRSTIMD = FIRGTILIOFIASGT (Md 169010359712. GOOATAIS) IN ID1410 
920 16 € 12297 > 33005603340 1891420 


CUMMAREAAAAAAAAARAAXARKAERARLARAAAKARARLA RURECAALAAAEAAAARAKAAKK AKARALAT 10 JO1930 


c 


TZZIT 15 NEGATIVE FOR STANDING MODES AND POSITIVE FOR SPINNING NODES 1NI01440 
CMMRMNARAAANAAARASRERARAAAWAAAARARAN ARAN RENAN ANNA AN AARAK EARNS RE REA RARER 18014 SO 
330 VT © VOTMATAI J) 1.901460 
CvT © COstvT) INJ01670 

Sv? = Simtvt) 14901480 

GO 10 330 19301490 

900 CVT © 100 1401300 
SVT © 100 14901510 
990 FP 2 100 14901920 
FR © 100 14 JB1$90 

FT © 100 183013840 

GO 10 (42003960)e OSCRE 14301930 
360 PO © ABSIPOOSS I YNECVT) 14301960 
VO © ABSIVOOSDSI VNECVT) 1m301970 

WO = ABSIVOOSS IyNeveSsvT /R) 1mg01980 

A © -30146159 14501990 
een 14301600 

17 6 THUR 37023800370 7991630 
$70 SAVEP © 1.0713 0161990P00TFLA) f.. M1420 
CALL INTER( Ae oPS! ot) 14JD1830 
CPST + COStPSSD 14016600 
19x @ POSCPS! 1NJD16$0 
CALL INTACIPPEI D> OPPELYs TOK. CPR 1%901660 

F © OPROCPS! 14581670 
COLL INTESIF y INTEG EL etme 19503680 

CP © SAVEPOINTES 19303690 

seo te ¢ THER 39004000990 1™J01700 
390 SAVER © 2.0/13.141990VO0TFLR) 1HJ03710 
CALL INTGR(AcB oPS1 oMUNBR) 14901720 
COST © COSTPST) 1*J03730 

39x © yOeCrs! tes01740 
CALL INTO CIPRELD© CPRELD> LPK> CPX) taso17se 

* « OPxeCPS! 1HJO1760 
CALL INTGS (Fe INTEGs El, rent INID1770 

FR = SAVERGINTEG 1mJ017600 
200 IF ( THAT» 41004200610 1"JO3790 
410 SAVET © LoO/(FoIQISSOVTOTFLT? 1MJB 1600 
CALL IMTER( A cB PSS oeUBR) 1302030 

CPs © COSIPS!) 14301620 

19m @ wOeCPSs 1WgD1830 
CALL INT CLPTER Ds OPTELI® fOxXe OPK) 1m 301840 

F © OPxECPSI 14301880 
CALL IMTGSEMSINTEGsE 1. 00m) . 14JB1860 
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COE ae wei 


INPUT ROUTINE: 
C1538 )--- - READ CHAMBER, NOZZLE DATA 
FOR NEW CASE 


CRUN = C RUN OR B RUN OR A RON 


7 
T 
AV HAMBER DATA ON SCRATCH TAPE 


(Uy) AS THE DATA AREA “ILL BE “RITTEN 
OVER BY A NES" LINK. 


: 


PLACE FREQUENCIES FOR DESCRIRTNG 
FUNCTION INTO UrrER COM‘ON 
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DECK = CHAMBER 
2 of 5 


RESTORE CHAMBER DATA ‘HICH 
FAS SAVED BEFORE ENTRY TO 
INJCTR PROGRAM. 


HOVE As Bays Cay FROM UPPER COMMON 
TO DATA BLOCK FOR "TRANS" OR "DDD" 


MOVE OUTPUT FROM "CCC" INTO MOVE OUTPUS FROM "CCC" 
INPUT BLOCK FOR "TRANS" INTO INPUT BLOCK FOR "LONGL" 
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CALCULATES HTR, HIT IF NOT H RUN; 
EXPANDS TABLE TO 0 POINTS 


CALCULATES (np, 7) FOR 
EACH FREQ. USING HT, HTT 
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LsCK = WGEN 


LONGITUDINAL MODE IF SNH IN 
PUT NEGATIVE; GENERATE FRE- 
QUENCIES ABOUT ABSOLUTE VALU 


AND SET SNH = 0.0 


GENERATE TSN FREQUENCIES 


O9S<W<sl.l S 


INDICATES FREQUENCIES 
W 2) 2910 OPS SS Seo HAVE BEEN GENERATED 
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i 
: 
; ENTRY = ‘TRANS 
DECK = HYMN 
4 lof 2 


CALCULATE LOW-ORDER TABLES 


Z, U (z), w, (2), (2) 5p, (z 


gCALCULATE HIGH-CRDER TABLES 


: ds He ir Ge a > ete. 
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2 of 2 (2° ) 


BEGIN ee eee EACH FREQUENCY: 
Suh - * 


eo" 


EVALUATE 6 INTEGRALS BY SIMPSON's 

(QR BOOLE'S) FORMULA; 

CALCULATE FIRST-ORDER SOLUTION Ho» H 
I 


+3 


CALCULATE SECOND-ORDER SOLUTION 
HR, HTT 
INCLUDING INJECTOR EFFECTS 


END OF FREQUENCY LOOP 


PRINT OUTPUT: 
Vy Hy H 
OR 


ne HL» Hs Hop » Ener 


=> 200 


DUMP TABULAR FUNCTIONS 


« 


Te ee ee 


F 


FOR EACH FREQUENCY, “’, 
CALCULATE HIR, HTI 


SAVE ONLY THOSE ENTRIES 
FOR WHICH HTR > 0 
(ERROR RETURN IF NONE) 


EXPAND TABLE 
(BY INTERPOLATION) 


TO hO POINTS 


<> 


>10 | 


PRINT OUTPUT TABLE | 
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VEIOCITY POTENTIAL DESIRED MACH MO. MACH MO. WILL BE 
TABLE AMD MACH MO. IS INPUT TO CALCULATED BY 
ARE BOTH INPUT °TBLCAL” *TBLCAL* 


SOLVE EIGHT SIMULTANEOUS 
DIFFERENTIAL EQUATIONS BY 
ADAMS METHOD (SUBROUTIXE 
"ADMINT") 


t 


DO FOR EACH COMPUTE AND STORE NOZZLE 
FREQUENCY ADMITTANCE COEFFICIENTS 


FOR EACH OF Kili STATTOHS 7, 
CALCULATE RADIUS R (2), 

ALD AREA A (Z) FROM THOUT 
iNZZLE GEOMETRY. 


CALCULATE TARLE A? VS, AM 
(AREA VS. MACH HNO.) 

FROM ISENTROPTC FLO’ 

RELA TTONS 


TOR FACH A (Z) THTERPOLATE TN 
A’, AM TABLE TO FTN) CORRES- 
PONDING MACH LUMBER. 


YETERNINE SQUARE OF REDUCED 
VELOCITY FROM MACH NUMBFR 


CALCULATE VELOCT™Y POTENTIAL ,XX 
= 2K z 


wa? 
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emRy = | LEGL 
Ycr = > 


DETEAAIME EFFECTIVE STA°TCu, AX, 
FO COiCAFPRATED CORBUSTTON 


(TIBERT FOUR SPEOQUTCIES 
REMEZ. PACH "HOUT “aT? ) 


ow tN = od 
DETR Oot, FOP BAH 


PREAH NOY : 


FOR ZACH FREQUENCY, CALCULATE 
AD RIT gs 7 1 


BUILD TABLE OF d #%& vs ~™ 
aw 


INTERPOLATE FOR 
ee (a2 = C) 


: dw 


: DETERMINE Hop (',), Hoy (,,) 


CALCULATE 72 (‘1 ), Tos ae 


PRINT MINIMUM VALUES 


<> 


READ TAPE 13 
(DATA FRO; PREVTGUS CASE 
HAS REci. SAVED Gii SCRATCH 
TAZ (13) 


Ce 138) este, READ INJECTOR 
DATA 


SAVE DATA OL SCRATCH TA™E 
(13) AS 2ZART OF INPU™ APFA 
AY BE DESTROYED 


ZERO STORAGE FO: 
SJECTO: TH.UT 


RETURN T) CHAIR?! 
VTA ATHY 


RETURN 


= 


AR ee Bp | 


RECOVER INPUT DATA FROM 
SCRATCH TAPE (13) 
CALCULATE “C’FOR EACH 
ELEMENT 


SAVE -“: FOR "INJDIS" 
(SCRATCH TAPE 13) 


2 DETERMINE —: (r) 


FOR EACH RADIAL BAND 
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BEGIN LOOP FOR EACH 
OF NE E.EMENTS 


1) 


CALCULATE DENOMINATOR D 
<0 STANDING <> > 0—SPINNING 


CVT=1.0 
SVT#1.0 


2 SOREN Ene es oe ee le ee 


DECK = INJDS 
2 of 2 


ACCUULATS TUIDRATORS PCR 
Ayers Pym Car 


RINTY | PF 
POR sg 
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END OF ELEMENTS 100 "aTIP: 
FOR NEXT ELITN™ OP ccnurtving 
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